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The surprisingly rapid relaxation of the sustainable current density in the critical state of single 
crystalline Bai_ x K ;E Fe2As2 is investigated for magnetic fields oriented parallel to the c-axis and to 
the afo-plane respectively. Due to the inadequacy of standard analysis procedures developed for flux 
creep in the high temperature superconducting cuprates, we develop a simple, straightforward data 
treatment technique that reveals the creep mechanism and the creep exponent fi. At low magnetic 
fields, below the second magnetization peak, varies only slightly as function of temperature and 
magnetic flux density B. From the data, we determine the temperature- and field dependence 
of the effective activation barrier for creep. At low temperatures, the measured current density 
approaches the zero-temperature critical current density (in the absence of creep) to within a factor 
2, thus lending credence to earlier conclusions drawn with respect to the pinning mechanism. The 
comparable values of the experimental screening current density and the zero-temperature critical 
current density reveals the limited usefulness of the widely used "interpolation formula". 

PACS numbers: 74.25.Ha,74.25.0p, 74.25.Wx 



I. INTRODUCTION 

Recently, the measurement of the critical current den- 
sity j c of superconducting iron-based compounds has 
been recognized as a useful tool for the characteriza- 
tion of microscopical and nanoscaleP^ disorder in these 
materials. However, the study of vortex pinning and 
the critical current density is compromised by surpris- 
ingly large thermally activated flux creep^^l m the Bean 
critical stateP^ In some materials, such as single crys- 
talline (Ba,K)Fe2As2, the logarithmic creep rate S = 
—dhuj/dhit of the sustainable current density j ap- 
proaches that previously measured in some of the high 
temperature cuprate superconductors PED Sizeable creep 
rates influence the magnitude, and, potentially, the tem- 
perature and flux density-dependence j(T,B), thereby 
compromising the analysis of fundamental vor tex pinning 
mechanisms in the iron-based superconductors^! and the 
understa nding o f possible phase transitions of the vortex 
ensemble.™! 

A detailed analysis of flux creep in the iron-based su- 
perconductors is therefore justified. In the cuprates, such 



analysis has unveiled the non-logarithmic nature of vor- 
tex creep, a direct consequence of the non-linearity of the 
relevant potential barrier U(j) opposing thermally acti- 
vation as function of the driving force Bj. In turn, this 
nonlinear behavior arises from the elasticity of the vortex 
ensemble, i.e., the fact that this can be deformed continu- 
ously on very different length scales. In contrast to single- 
particle creepp^l the relevant activation barrier does not 
depend algebraically on the driving force, but, rather, in- 
creases steeply at low driving forces because of the nonlin- 
ear increase of the size of the critical nucleus EMUI i n g en _ 
eral, the size of the critical nucleus increases as an inverse 
power-law in j, leading to the well-known relatio n 15 " 17 ! 

u(j)=u c (^y. a) 

The value of the creep exponent /i depends on the di- 
mensionality of the critical nucleus as well as that of the 
elastic manifold as a wholeP^ In the case of single crys- 
talline iron-based superconductors to be considered here, 
the latter can be either a single one-dimensional (ID) 
vortex line, or the 3D vortex ensemble. The law (TTh does 



2 



not extrapolate to a zero activation barrier at large driv- 
ing force; one therefor e freq uently resorts to the so-called 
interpolation formula! 16 * 17 * 



U(j) = U c 



(2) 



Among the various me thods to experimentally estab- 
lish the £/(j)-dependencej 10 l 18 * ^the analysis of the loga- 
rithmic time dependence of 1/5 is the most reliablePH22l 
for Eq. ([2]). one has 1/S oc /iln[(t + t)/r], with t Q a 
constant determined by transients at the onset of relax- 
ation, and r = (Aj c /Eo) {ksT /U c ) a normalization time 
determined by the barrier magnitude, a factor E$ (with 
dimension of electric field) related to the details of the 
creep mechanism, and a sample "in ducta nce" A = hqo? 
(a is a relevant sample dimension) ) 23 l 24 l An alternative 
"dynamic" method based on the dependence of the sus- 
tainable current density on the sweep-rate H a of the ap- 
plied magnetic field in magnetic hysteresis loop measure- 
ments was used in Refs. OB 20, and l2"5Tf2"Tl In princi- 
ple, the product ( T/Q)d \nj/dT, with Q = d\nj/d\nH a , 
directly yields ^l 5 l 6 l 25 tf Hl The potential of local mag- 
netic measurements of the mag netic induction B was 
exploited by Abulafia et al.pSl who reconstructed the 
current-voltage characteristics of YBa 2 Cu3C>7_,5 single 
crystals from Ampere's law j t 
rent density and Faraday's law 



(U 1 dB jdx for the cur- 



E(x) 



{dB/dt)dx' 



(3) 



for the electric field (integrating from the sample centre 
to its perimeter) . 

As opposed to the cuprates, the dynamic range 
over which creep data can be collected in the iron- 
based superconductors is insufficient to reliably deter- 
mine S(lnt), Q(\nH a ), or the curvature of E(j) = 
E exp(—U(j)/kBT). Moreover, the factor d\nj/dT oc 
[1 — (ksT /U c ){dU c / dT)] used in the dynamic method ad- 
mixes the temperature dependence U C (T) with the creep 
exponent ji. As a result, it is difficult to distinguish be- 
tween a logarithmic U ~ U c ln(j c /j) with a temperature- 
dependent U c , and a so-ca lled "negative //"-type barrier, 
U = U c [l - (j/j c )\»\]WE®One should therefore resort to 
ot her const ructions, such as that proposed by Maley et 
al. j 5 l 13 l 19 ]Hl which combine the results of relaxation mea- 
surements at different T. 

Below we show that even though this method is fraught 
with shortcomings, it can be suitably adapted to yield 
reliable results. In particular, a plot of the average 
(— ksTdlnQdj /dt\)/dj) vs. the average (j) unambigu- 
ously yields the curvature of the £/(j)-relation and there- 
fore the flux creep mechanism. Applying this to single 
crystalline (Ba,K)Fe2As2, one finds, first of all, that the 
creep rate is insufficient to qualitatively modify the field 
dependence j(B) as obtained from magnetic hysteresis. 
Second, the low-temperature sustainable current density 
j is, typically, less than a factor of 2 lower than the value 
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FIG. 1: (Color online) (a,b) In-phase fundamental (T' H ) and 
amplitude of the third harmonic ac transmittivity ( | Tj^3 1 ) as 
function of the temperature, at an applied field (j,oH a = 1.5 T, 
for different indicated frequencies of the ac field, (c) Irre- 
versibility field Hi rr , as determined from the vanishing of 
\Tzh\, indicated by the arrow in (b), for H a \\ c and H a \\ ab 
respectively. 



the critical current density j c would have in the absence 
of creep. Finally, the obtained /i-values indicate that the 
same creep mechanism is relevant for all temperatures 
up to the transformation of the vortex ensemble at the 
so-called "second magnetization peak".E2l We discuss the 
/i-value and its increase as function of magnetic field in 
terms of the interplay of strong pinning by nm-scale het- 
erogeneities and weak collective pinning by atomic scale 
point defects in the material. 



II. EXPERIMENTAL DETAILS 

The Bao.6Ko.4Fe2As2 single crystals with T c ~ 38.1 
K (Fig. [lj were grown by the self-flux method, using 
FeAs as the self-flux. Details of the growth can be found 
in Refs. [32 33. For the experiment, crystals were cut 
to regular rectangles of dimensions 350 x 160 x 70 /im 3 
using a wire saw. The spatial distribution of the lo- 
cal induction perpendicular to the surface of the crys- 
tals in the critical state was measured using an array 
of microscopic Hall sensors, fashioned in a pseudomor- 
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FIG. 2: (Color online) (a) (o) Hysteresis loop of the gradient 
dB/dx oc j° b versus the local induction B, measured perpen- 
dicular to the long edge of the Bao.6Ko.4Fe2As2 single crystal, 
for H a || c at T = 10 K. The sweep-rate of the applied mag- 
netic field was 0.8 mT/s for the virgin magnetization loop (up 
to fJ,oH a = 0.5 T), while further measurements were carried 
out with a sweep-rate of 4 mT/s. The four measurement se- 
quences ( □, O, A, •) illustrate the decay of the flux gradient 
at different applied fields, over a period spanning 5 to 5000 
s after field-cooling, and reduction of the external field H a 
to the measurement field H a — AH a (with AH a = 0.2 T). 
(b). Time dependence of the local induction gradient for an 
applied field of 0.6 T and T = 8 K. 



phic GaAlAs/GaAs heterostructure using ion implanta- 
tion. The 10 Hall sensors of the array, spaced by 20 /im, 
had an active area of 3 x 3 /im 2 , while an 11 th sensor 
was used for the measurement of the applied field. For 
measurements with field parallel to ab, an array of sen- 
sors spaced by 10 /im was used. The Hall sensor array 
was placed on the center of the crystal surface, perpen- 
dicular to the long crystal edge, and spanning the crystal 
boundary. In this manner, hysteretic loops of the spa- 
tially resolved local induction were measured as function 
of the applied magnetic field H a , and as function of tem- 
perature. The sweep-rate of the applied magnetic field 
was 40 G/s (noH a = 4 mT/s). In all experiments, the 
applied magnetic field was not only much smaller than 
the upper critical field H C 2, but also smaller than the 
field at which the second magnetization peak occurs.^ 

Under all circumstances, the profiles of the flux density 
B were well described by the critical state modelj^ see 
Fig. [3] allowing for the straightforward extraction of the 
local screening current density j ps (2/ HojdB jdx. The 
factor 2 takes the finite dimensions of the crystal into 
account; namely, the field gradient on the end surface of 
a semi-infinite bar in the critical state is half the field 
gradient in the interior.^ Measurements were done with 
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FIG. 3: (Color online) (a) Relaxation of the critical state 
flux profile in the Bao.6Ko.4Fe2As2-crystal, for H a || c, at 
T — 8 K. Open symbols show the flux density profile across 
the Bao.6Ko.4Fe2As2 single crystal for four successive times, 
while closed symbols illustrate the flux density profiles ob- 
tained after the annealing of the critical state to the indicated 
temperatures, (b). Flux-creep activation barrier as function 
of the screening current density, as extracted using Eq. @. 



the Hall array parallel to the aft-plane and H a | c, a 
configuration that measures the usual screening current 
density in the aft-plane, j — j c ah , and with the array || c 
for H a || ab, which yields the aft-plane screening current 
density corresponding to vortices moving parallel to 
the c-axis (perpendicular to the FeAs layers). 

Measurements of the ac screening were performed with 
the same set-up. A sinusoidally time- varying field of mag- 
nitude 1 Oe and frequency / is applied colinearly with 
the dc field. The ac component B ac (f,T) of the local 
induction is then measured using the Hall probe array. 
Results, presented as the in-phase fundamental ac trans- 
mittivity T' H = [B ac (f, T)-B ac (f, T « T c )]/[B ac (f, T » 
T c ) - B ac (f,T < T c )\ and the third harmonic \T 3H \ = 
B ac (3f,T)/[B ac {f,T » T c ) - B ac (f,T < T c )]^ are 
shown in Figs, [T^t, and [TJ) respectively. The non-zero 
value of the latter signals the existence of a non-zero 
critical current density j c ; the vanishing of B ac (3f) at 
high-temperature is used to trace the dc irreversibility 
field Hi rr (T), for dc field aligned along the c^axis and 
the a6-plane respectively (Fig. [l];) . 

As for the magnetic relaxation experiments, these were 
carried out on the decreasing field branch magnetization 
branch only (corresponding to flux exit relaxations), in 
order to prevent possible influence of surface barrier re- 
laxation. The external magnetic field was applied at a 
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temperature T > T c , at which the Hall probe array was 
calibrated (with respect to the applied field) . The crystal 
was subsequently field-cooled to the measurement tem- 
perature T e < T c , and the field reduced by an amount 
AH a to the measurement field. Care was taken that AH a 
exceeded the the field of full flux penetration, so that a 
full critical state is established. After waiting 5 s to allow 
for the settling of the magnet, the flux density values at 
the different Hall sensor positions were measured as func- 
tion of time, for a period of 5000 s. Furthermore, flux 
creep annealing experiments^ were performed by heat- 
ing the sample to T a > T e , returning to the experimental 
temperature T e , and re-measuring S over a period of 600 
s. This procedure is equivalent to performing measure- 
ments at the effective time t = t exp[U(j(T a ))/k B T e ], of 
the order 10 6 s. 



III. RESULTS 

Fig. [2^ a) shows hysteresis loops of the flux density 
gradient dB/dx as function of the local induction B 
for H a || c, at T = 10 K. The width of the loop is 
proportional to the screening current density j c ah in the 
a6-plane, for field || c. It has the characteristic shape 
found in all charge-doped iron-based superconductors.^ 
a "central peak" of j around B = 0, followed by a drop 
j cx i?" 1 / 2 characteristic of strong vortex pinning by 
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FIG. 4: (Color online) Treatment of creep data for fj,oH a — 0.2 
T || c, using the method of Ref. 1291 The data at the bot- 
tom of the frame show the ^'-dependence of the product of 
T e and the logarithm of the electric field, such as extracted 
using Eq. j3t, for different experimental temperatures T e (left- 
hand axis) . The master curve at the top (righthand axis) was 
compiled from these data using a temperature- and current 
density-independent Eq; however, as in Ref. 1291 Eq was as- 
sumed to be proportional to B; this yields a /x-value of 0.8. 



nm-scale point defects^* 34 * 35 ! or heterogeneities. As B 
increases, the strong pinning-contribution to the criti- 
cal current becomes irrelevant, and j(B) saturates to a 
field-independent value determined by weak collective 
pinning of individual vortices! 3 ^ presumably by the K 
dopant atoms in the material!^ The magnitude of j c ah is 
very similar t o that f ound in other measurements on the 
same material ! 1 * 38 * 39 ! 

At all magnetic fields, the sustainable current relaxes 
as function of time, with S ~ —0.06 [Figj2^b)]. Creep is 
logarithmic in time, with curvature indicative of a non- 
linear U(j) relation. However, contrary to the cuprate 
superconductors, and as illustrated by the near-linear 
evolution of the experimentally determined flux-creep ac- 
tivation barrier in the inset to Fig. [3] the accessed dy- 
namical range of S(j(t)) is too small in the iron-based 
superconductors to reliably extract a value of \i directly. 
The effect of flux creep-annealing is also illustrated in 
Fig. [3] which shows the relaxation of the trapped flux- 
profile at T — 8 K together with profiles obtained after 
annealing at various temperatures. 

The time-dependent current density is analyzed using 
the method of Maley et aZ.P^This yields the experimen- 
tal current density-dependent activation barrier as 



Ue(j) 



-k B T e \n 



(4) 



with c a temperature-independent constant, to be cho- 
sen so that segments corresponding to the U e (j) relation 
measured at different temperatures line up to form the 
U(j) relation representative of the vortex creep mecha- 
nism governing the magnetic relaxation. For the collec- 
tive creep mechanism of Ref. [TBI 



Jc 



{(k B T/U c )ln[(t + t) /t]}Vm' 



so that the method giveiP^ 



k B T e In 



U(j)-k B T e ln 



k B T e j 
U(j) [IT 



(5) 



(6) 



hence, c = k B In [k B T e j/U(j)fiT] actually depends log- 
arithmically on temperature. For a put ative logarith- 
mic dependence U — U c ln(j' e /j) j 4Q l 41 l yielding j — 



j c [(t Q + t) /t] 



fc/sT fi In 



-k B T/U c 



haiP 



= U(j) - k B T e In 



k B T e j c 
U e T 



(7) 



The following shortcomings of the method of Ref. [TH] 
are identified. While, at low temperature, only a sin- 
gle c-value will satisfy the requirement of lining up the 
measured segments, at higher T the temperature depen- 
dence of c becomes important. Second, the prefactor 
U C {T) ee U c (0)f(T), as introduced in Eqs. and 
itself introduces a more important temperature depen- 
dence. This can be recovered by dividing the result of 
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FIG. 5: (Color online) Compiled double-logarithmic plot of 
(d[—kBT e dln(\dj/dt\)]/dj) vs. (j), as suggested from Eq. |9|, 
for all temperatures, and the two orientations (|| c and || afc) 
of the applied magnetic field. 



Eq. Q at each T e by phenomenological factors f(T e ), 
with, again, a degree of arbitrariness. Correcting for the 
temperature dependences of c and U c in different ways re- 
sults in different final results for the compiled U (j)— curve. 
Finally, the method supposes that the same U(j) mech- 
anism governs flux creep at all the temperatures used to 
reconstitute the experimental U (j)— curve. In the present 
set of experiments, we find that applying different pro- 
cedures to cope with these different temperature depen- 
dences yield an error bar on fx that amounts to 50 to 
100 % of its value. 

Because it involves the integral |3| rather than the lo- 
cal value dj/dt, the method of Ref. 1291 has the merit of 
yielding more accurate [/(j)-values. In our experiments, 
we determine the electric field by integrating to the sam- 
ple boundary; in accordance, the relevant value of the 
current density is that which corresponds to the slope 
dB/dx at the boundary. Fig. [4] shows —k B ThiE thus 
determined, for various temperatures, and noH a = 0.2 
T || c. However, if one wishes to extract flux-creep ac- 
tivation energies, one is faced with the same arbitrari- 
ness concerning the factor c in Ref. W\ now contained 
by the factor Eq. Compiling a full [/(j)-curve using a 
temperature- and current-density independent Eq yields 
the illustrated master curve, suggesting that U(j) cx 
with fj, = 0.8. However, if one takes Eq to be propor- 
tional to Bj, as in Ref. [551 one obtains a curve that fits 
U(j)ocj-°- 25 . 

To avoid ambiguities, we evaluate the averaged current 
density-derivative (— k B T e d\n \E\/dj) at each T e . This 
procedure has the advantage of eliminating the prefactor 
c (or, equivalently, Eq) from the analysis. For Eq. (|T|), 
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FIG. 6: (Color online) Experimental values of the creep expo- 
nent [i for the two orientations of the applied magnetic field. 



one has 

.k B Td\u\E\ 



d.i 



= (fi+l)k B T e dU(j) 



A* — 

3c 



(J)J 



dj ' 
(U c » k B T) 



while for power-law cre ep^li3 



k B Td\n\E\ 
do 



Uc 

or 



(U c » k B T) (9) 



A double-logarithmic plot (Fig. |5| of the average 
(—k B T e dln\E\/dj) vs. (j) for all different temperatures 
T e therefore directly reveals the creep mechanism, as well 
as the value of fj,, from the deviation of the linear slope 
from 1. Whether the condition U c 3> k B T is verified 
can be easily checked from the Bean (straight-line) -like 
nature of the flux profiles in Fig. [3J 37 

From the slopes in Fig. [5] we find that Eqs. ([IJ and 
([2]) describe the data satisfactorily. The corresponding 
//-values are rendered in Fig. [6] For H a || c, fi = fi c » 
0.65 at low fields, slowly increasing towards /i c w 0.8 
at higher fields, while for H a i, \\ ab, [i ah w 1.5 at low 
fields, decreasing towards /i ob 1 at higher fields. The 
results for H a \\ c are remarkably similar to those found 
in Ba(Fe .92Coo.o8)2As 2 by Shen et al.W 

The obtained /i-values are checked by plotting 
— k B T e \n\E\ vs. for all temperatures and fields. 

Fig. [7] shows that this yields straight lines for both 
field orientations, as required. Deviations from linear- 
ity are only apparent below 8 K, and above 24 K, which 
means that /1 is T-independent for the greater part 
of the investigated temperature range. That is, the 
same mechanism governs vortex creep for all T < 24 
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FIG. 7: (Color online) Current-density dependence of the 
barrier U(j) for thermally activated vortex creep in the 
Bao.6Ko.4Fe2As2 crystal, for (a,b) fj,oH a = 0.6 T || c, fx = 0.75; 
and (c,d) ^ H a = 0.5 T || ab, n = 1.13. (b) and (d) 
show the (logarithm of) the raw electric field data, multi- 
plied by the experimental temperature, — fcsTedhi \E\, ver- 
sus (dB/dx)~^ oc (a,c) show the full, compiled curves 
of the creep barrier, obtained by adding the relevant factor 
fc_eT e lni5o for each temperature. 



K. The slopes in Fig. [7] yield U c j£, while the inter- 
cept with the abscissa corresponds to where j e = 
j c {U c /kBT e \n[k B T e j/U{j) l XT]} 1 '' 1 . With the uncer- 
tainty regarding fi removed, the curves obtained for dif- 
ferent T e can now be compiled into a "zero-temperature" 
curve by adding factors fcsTelnii'o, and dividing by 
/(T e ). For the lowest two temperatures, we set /(T e ) = 1 
to obtain an unequivocal Eq = 0.019. For higher T e , 
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FIG. 8: (Color online) Comparison of the temperature de- 
pendence of the activation barrier U c (T)[j c (T) / 'j c (0)] M (open 
symbols) extracted from the creep experiments, with that of 
the measured screening current density j(T) at the onset of 
relaxation (closed symbols). The j(T) values (closed sym- 
bols) are tantamount to those one would measure in magnetic 
hysteresis experiments. Large closed symbols on the lefthand 
abscissa denote the j c (T = 0)-values, determined from the ex- 
trapolation of the activation barrier to zero in Fig. [9] Drawn 
lines denote Eq. q5b , evaluated using the experimental values 
of the product Ucjg, and ln(to + t/r) = 21. Das 
guides to the eye. 



Dashed lines are 



the value of f{T e ) is adapted in order to line up the 
relevant data segments to obtain a continuous compiled 
curve, with continuous derivative. The result is given in 
Fig. [9] which shows the extracted activation barrier U(j) 
as function of j~ (top panel) and as function of j (bot- 
tom panel), for various magnetic fields || c. The extracted 
temperature dependence U c jj? for H a \\ c is depicted in 
Fig-i 



IV. DISCUSSION 

Inspection of the creep barrier in Fig. [9] shows devia- 
tions from power-law behavior at low temperature < 8 
K, and above 24 K. In all cases, the dependence U(j) at 
the lowest temperature appears to be linear in the current 
density j . This behavior is interpreted as being due to the 
proximity of the measured screening current density to 
the (pinning) critical current density. A linear extrapola- 
tion of the U(j) curves to zero (dotted lines in Fig. [9]) thus 
provides an estimate of the critical current density j c (0) 
in the limit T — > 0. The obtained values range between 
j c (0;0.6T) = 6.7xl0 9 Am" 2 toj c (0;H a = 0) = 9.7 xlO 9 
Am" 2 for field || c, and j c (0; 0.5 T) 2.3 x 10 10 Am~ 2 for 
H a || ab (Fig. [8]). These values can be factored out from 
the slopes of the curves in Fig. |9ja,b) to yield creep bar- 
rier values of U c (0)/kB = 60 - 85 K, decreasing with in- 
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FIG. 9: (Color online) (a) Plot of -k B T e d\n\E\ vs. 
(dB/d:r)~ M oc j~ M , for different magnetic fields strengths || c. 
(b) Flux-creep activation barrier compiled from relaxation ex- 
periments at different temperatures, for the same fields as (a). 
The dotted lines show the linear extrapolation of the low- 
temperature barrier to Jc(0). The drawn lines show fits to 
the collective creep theory, Eq. |l]). (c) Data for applied field 
HoH a = 0.2 T, showing the inadequacy of the interpolation 
formula |2| (dashed-dotted line). 

creasing values of field II c, see Fig.[8j and U c {0)/k B = 24 
K for H a || ab. The low-temperature values of the mea- 
sured screening current density j closely approach j c (0). 
These values therefore yield a good order-of-magnitude 
estimate of the bulk pinning force. 

The Inset, Fig. [9jc), compares the high- and low cur- 
rent limiting behavior of the current density-dependent 
creep barrier for /j,oH a = 0.2 T || c to the often-used 



interpolation formula, Eq. It is clear that the cur- 
rent density range over which the experimental barrier 
crosses over from the high-current, linear-in-j, to the 
low current-behavior given by Eq. is much narrower 
than what can be described using Eq. 0. In fact, the 
dashed-dotted curve shows that applying the interpola- 
tion formula ^ yields a gross overestimate of the pinning 
critical current density, as well as a possible overestimate 
of fx. It therefore seems imperative to use Eq. 0, which 
was derived on physical grounds; 1 ^ rather than the phe- 
nomenological formula ([2j. 

The extracted parameter values can be used to cross- 
check the analysis. The drawn lines in Fig. [8]render a nu- 
merical evaluation of Eq. [5] using the experimentally ex- 
tracted Udc products, and ln(to+t / t) = 21. In the inter- 
mediate temperature range, at which creep is described 
by the barrier 0, the agreement with the temperature- 
dependent screening current density values j(T) at the 
onset of the relaxation (such as these might be mea- 
sured during a field-sweep measurement) is more than 
satisfactory. Given that all times the current density is, 
to good approx imation , given by the equation U(j) — 
fc B rin(to + VT). | 23 ' 24 ' 41 | one has ln(i +t/r) = U(j)/k B T, 
where the numerator and denominator can be simply 
read from Figs. [4] and [9] At the onset of relaxation (t <~ 5 
s), one finds values ranging from 20.5 (for 4.2 to 12 K) 
to 52 (at T — 26 K). Thus, the value ln(t + t/r) = 21 is 
reasonable in the intermediate temperature range. The 
expression r = (Aj c / Eo)(kBT/U c ) ~ 10 -8 s allows one 
to estimate ln(io + f/r) independently; for i ~ 5 s one 
again has a value of 20. For measurements performed 
with a commercial superconducting quantum interfer- 
ence device-based magnetometer, t ~ 100 s, so that 
ln(io + t/r) is slightly larger. 

Fig. [8] shows that the field dependence of the screen- 
ing current density, as expressed by Eq.([5]), is contained 
by the parameter U c j£. Hence, the j (independence is 
not the consequence of a field-dependence of the creep 
process, which would be reflected by a strongly field- 
dependent /i; rather, it reflects the intrinsic field de- 
pendence of the pinning force. The _B -1 / 2 -dependence 
of the screening current density in iron-based super- 
conduc tors w as recently interpreted in terms of strong 
pinnin^2il25| nm -sized heterogeneities JD® l n this re- 
spect, the creep exponent /i ~ 0.6 -0.8 at low fields 
parallel to the c-axis, comparable to fi = 0.5 found for 
single vortex creep in the Bragg-glass phase in single 
crystals of the cuprate high temperature superconductor 
Bi2Sr 2 CaCu208 + 5j^2142l comes as somewhat of a surprise. 
Namely, vortex pinning at low fields in E^S^CaC^Os+a 
is thought to be not in the str ong pinning limit, but in the 
opposite, weak pinning limit !^ In the field- values under 
scrutiny, the field-dependence observed in Fig. |2]suggests 
that one is dealing with strong pinning (see also Ref. [T3J). 
A possible explanation for the similarity of the creep ex- 
ponents in the two cases is that the velocity of the flux 
lines at low and intermediate currents is limited by the 
progression of roughened vortex segments spanning the 
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region between strong pins through the background of 
weak pins. This would seem natural given that weak 
pinning in the iron-pnictide superconductors is thought 
to be due to the local fluctuations of the dopant atom 
density, while strong pinning would originate from inho- 
mogeneity of the dopant atom density on a much larger 
length scale. Thus, vortex segments would progress more 
or less continuously through areas of homogeneous dop- 
ing, before becoming stuck on a "strong pin" , that is, a 
region in which the dopant atom density averaged over 
several dozen nm significantly deviates from the overall 
mean. The critical current density, and the screening 
current density at low temperature, would then be de- 
termined by the strong pins, while creep at intermediate 
and high temperatures would be determined by the weak 
pinning background. An alternative hypothesis is that 
the creep exponent would be determined by the shape of 
the energy distribution of the vortices pinned by large- 
scale heterogeneity, much in the same way as this was 
proposed for creep through columnar defects in the so- 
called variable-range hopping regime.^ 

Finally, we remark that the consistently larger \i- 
values found for H a \\ ab are not unexpected, because, 
in the investigated orientation, j \\ ab but _L H a . Vor- 
tex lines are therefore oriented in the a6-plane, but are 
forced to move parallel to the crystalline c-axis. This is 
the hard direction for vortex motion, requiring the nucle- 
ation of vortex loops in the c-direction, a process limited 
by the value of the out-of-plane vortex lattice tilt modu- 
lus c^~ 4 ~ £^" 3 C44. This exceeds the non-local tilt modulus 
£44 for vortices 1 1 c by the inverse cube of the anisotropy 
factor e\ ~ 0.4. 

V. SUMMARY AND CONCLUSIONS 

Iron based-superconductors show a surprisingly large 
value of the relaxation rate S of the irreversible magneti- 
zation. Since the relaxation rate is too small to directly 
extract the current-density dependence of the flux-creep 
activation barrier from the time dependence of l/S", an 



alternative, straightforward analysis method is proposed, 
that directly yields the relevant creep mechanism, and 
the temperature dependence of the pinning parameters in 
the absence of creep. Applying this to single crystalline 
Bao.6Ko.4Fe2As2, we find evidence for nucleation-type 
(collective) flux creep, with a weakly field-dependent 
creep exponent, fi = 0.6-0.8 for magnetic fields oriented 
along the c-axis, and a slightly larger /i = 1.2 - 1.5 for 
field along ab. Several hypotheses leading to such /1- 
values, among which, the combined action of strong and 
weak pinning centers, and a non-trivial pinning energy 
distribution function are proposed. At low temperature, 
the screening current approaches the pinning critical cur- 
rent; therefore, meaningful information on flux pinning 
in the iron-pnictide superconductors can be directly ex- 
tracted from low-temperature (T < 5 K) magnetic hys- 
teresis experiments. The field dependence of the screen- 
ing current density is found to arise from the underlying 
mechanism of pinning, and not from varying creep rates 
due to the flux creep process. Finally, it is found that the 
crossover between low-current and high-current behav- 
ior of the flux creep activation barrier is poorly described 
by the so-called "interpolation formula". 
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